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g : Abstract 
^ ■ 

The abihty of the Rigged Hilbert Space formahsm to deal with continuous 
^ • spectrum is demonstrated within the example of the square barrier poten- 

tial. The non-square integrable solutions of the time-independent Schrodinger 
equation are used to define Dirac kets, which are (generalized) eigenvectors of 



^ , the Hamiltonian. These Dirac kets are antilinear functionals over the space of 

physical wave functions. They are also basis vectors that expand any physical 
wave function in a Dirac basis vector expansion. It is shown that an accept- 
^ ' able physical wave function must fulfill stronger conditions than just square 

. integrability — the space of physical wave functions is not the whole Hilbert 

space but rather a dense subspace of the Hilbert space. We construct the 
position and energy representations of the Rigged Hilbert Space generated by 
the square barrier potential Hamiltonian. We shall also construct the unitary 
operator that transforms from the position into the energy representation. 
We shall see that in the energy representation the Dirac kets act as the anti- 
linear Schwartz delta functional. In constructing the Rigged Hilbert Space of 
the square barrier potential, we will find a systematic procedure to construct 
the Rigged Hilbert Space of a large class of spherically symmetric potentials. 
The example of the square barrier potential will also make apparent that the 
natural framework for the solutions of a Schrodinger operator with continuous 
spectrum is the Rigged Hilbert Space rather than just the Hilbert space. 



Typeset using REVT^ 
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1. INTRODUCTION 



In the late 1920s, Dirac introduced a new mathematical model of Quantum Mechan- 
ics based upon a uniquely smooth and elegant abstract algebra of linear operators defined 
on an infinite dimensional complex vector space equipped with an inner product norm [|l|. 
Dirac's abstract algebraic model of bras and kets (from the bracket notation for inner prod- 
uct) proved to be of great heuristic value in the ensuing years, especially in dealing with 
Hamiltonians whose spectrum is continuous. 

The Hilbert space (HS) was the first mathematical idealization proposed for Quantum 
Mechanics 0. However, as von Neumann explains in the introduction to his book 0], the 
HS theory and Dirac's formalism are two different things. Although there were attempts 
to realize the Dirac formalism in Hilbert space, there were a number of serious problems 
resulting from the fact that the Hilbert space cannot allocate such things as bras, kets, the 
Dirac delta function or the Dirac basis vector expansion, all of which are essential in any 
physical formulation of Quantum Mechanics that deals with continuous spectrum. Indeed 
in his textual presentation Dirac himself states that "the bra and ket vectors that we 
now use form a more general space than a Hilbert Space" (see page 40). 

In the late 1940s, L. Schwartz gave a precise meaning to the Dirac delta function as a 
functional over a space of test functions (cf. 0). This led to the development of a new 
branch of functional analysis, the theory of distributions 0. About the same time, von 
Neumann published the theory of direct integral decompositions of a Hilbert space induced 
by a self-adjoint operator (also valid for more general cases). This spectral theory was 
closer to classical Fourier analysis and represented an improvement over the earlier von 
Neumann's spectral theory 0. 

I. Gelfand always thought that von Neumann's spectral theory was not the whole story 
of the theory of linear operators defined on infinite dimensional vector spaces. Prompted by 
the theory of distributions, he and his school introduced the Rigged Hilbert Space (RHS). 
Starting out with this RHS and von Neumann's direct integral decomposition, they were 
able to prove the Nuclear Spectral theorem (also called Gelfand-Maurin theorem) [Q. This 
theorem justifies Dirac basis vector expansion. 

One of the aspects of Dirac's formalism, the continuity of the elements of the algebra of 
observables, was discussed in the early 1960s in Refs. If two operators of the algebra 

of observables satisfy the canonical (Heisenberg) commutation relation, at least one of them 
cannot be continuous (bounded) with respect to the Hilbert space topology. In Refs. [^,0, it 
is shown that there are subdomains of the Hilbert space that can be endowed with (locally 
convex) topologies that make those operators continuous; the largest of those subdomains 
is the Schwartz space. 

In the mid 1960s, some physicists []8|-|lO|] independently realized that the RHS provides 



a rigorous mathematical rephrasing of all the aspects of Dirac's formalism. In particular, 
the Nuclear Spectral theorem restates Dirac basis vector expansion along with the Dirac 
bras and kets within a mathematical theory. Later on, the RHS was used to accommodate 
resonance states (Gamow vectors) (cf. |]1T]-|T^ and references therein). Applications of the 
RHS formalism to Quantum Mechanics can now be found in some textbooks |[T8|JT9(| . 

Although there are some explicit examples of RHS in the literature (see for instance [pO[|), 
no example of the RHS generated by a Schrodinger Hamiltonian with continuous spectrum 
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has been constructed yet. Here we try to fill in this gap [|T7 . 

The dynamical equation that governs the behavior of a quantum system at any time is 
the time-dependent Schrodinger equation, 

th^^^it) = H^{t) , (1.1) 

where H denotes the Hamiltonian of the system and (p{t) denotes the value of the wave 



function (p at time t. In order to solve (|1.1| ), we associate to each energy E in the spectrum 



Sp{H) of the Hamiltonian a ket \E) which is an eigenvector of H, 

H\E) = E\E), EeSp{H). (1.2) 
These eigenkets form a complete basis system that expands any wave function ip as 

ip= [ dE\E){E\ip)= [ dEip{E)\E). (1.3) 



The time- dependent solution of Eq. ( |1.1|) is obtained by Fourier-transforming the time- 
independent solution (|1.3|) , 

(^(t) = j dE e~*^*/^ p}{E) . (1.4) 

If the spectrum of the Hamiltonian has a continuous part, and if the energy E belongs 
to this continuous part of the spectrum, then the corresponding eigenket \E) that solves 
Eq. ( |1.2| ) is not square integrable, i.e., \E) is not an element of the Hilbert space. Therefore, 
the eigenket \E) cannot represent an experimentally preparable physical state. 

We will show that the expansion ( |1.3| ) is only valid for those if that belong to a space of 
test functions ^ G 7i. We will also show that the kets \E) can be understood mathematically 
as continuous antilinear functionals over the space of test functions i.e., I^') G 
According to the RHS mathematics, equation ( [1.2|) means that 



{Hip\E) = E{ip\E) , for every ^ G $ , (1.5) 

where H is (the restriction to $ of) the self-adjoint Hamiltonian operator, which is a con- 
tinuous operator on the linear topological space For every such an operator, one defines 
the conjugate operator on by 

{Hcp\F) = {ip\H''F) , for all G F G (1.6) 

The operator if ^ is a uniquely defined extension of the Hilbert space adjoint operator 
(which for the case of a essentially self-adjoint operator coincides with the closure of H). 
Using the definition ( |1.6|) , we write (|1.5|) formally as 

H''\E) = E\E), |F)G$^ (1.7) 

which is understood as a functional equation over the space The quantities E and \E) 
are called generalized eigenvalues and generalized eigenvectors, respectively. 
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The general statement of the Nuclear Spectral theorem just assures the existence of the 
generalized eigenvectors \E), but it does not provide a prescription to construct them. In 
this paper, we construct the generalized eigenvectors \E) of the square barrier Hamiltonian 
along with the RHS. We shall use the Sturm-Liouville theory (Weyl theory) |^ to find the 
RHS of the square barrier potential. 

By applying the Sturm-Liouville theory to the Schrodinger equation of the square barrier 
potential, we will obtain a domain 'D{H) on which the Hamiltonian is self-adjoint. The Green 
functions, the spectrum, and the unitary transformation that diagonalizes our Hamiltonian 
will be also computed. The diagonalization of the Hamiltonian will allow us to obtain the 
energy (spectral) representation and the direct integral decomposition of the Hilbert space 
induced by our Hamiltonian. We will see why this direct integral decomposition is not 
enough for the purposes of Quantum Mechanics and why the RHS is necessary. Next, we 
will construct the space $ and therewith the RHS of the square barrier potential, 

$C7^C$^. (1.8) 

Dirac kets will be accommodated as elements of and the Schwartz delta function will 
appear in the energy (spectral) representation of the triplet ( |1.8| ). The Nuclear Spectral 
theorem will be proved, and it will be shown that this theorem is just a restatement of the 
heuristic Dirac basis vector expansion. 



2. STURM-LIOUVILLE THEORY APPLIED TO THE SQUARE BARRIER 

POTENTIAL 

2.1. Schrodinger Equation in the Position Representation 

In order to calculate the set of real generalized eigenvalues of the square barrier Hamil- 
tonian (the physical spectrum) and their corresponding generalized eigenvectors, we solve 
equation ( |1.7] ) in the position representation. 



{x\H'^\E) = E{x\E) . (2.1) 
The expression of the Hamiltonian in the position representation is 

{x\H'" \E) = f + V{x)] {x\E) , (2.2) 



^ 2m 

where A is the three-dimensional Laplacian and 



< r < a 

V{x) = V{r) ={Voa<r<b (2.3) 
b < r < oo 



is the square barrier potential. Writing Eqs. (p.l|) and ( p.2|) in spherical coordinates and 
restricting ourselves to the case of zero angular momentum, we obtain the radial time- 
independent Schrodinger equation. 
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+ V{r) ) x(r; E) = Ex{r; E) . (2.4) 



2m (ir^ 

Thus our Hamiltonian in the radial representation is given by the differential operator 

h^-^^,+V{r). (2.5) 

Throughout this paper, the symbol h will be used to denote the formal differential operator 
( p.5| ). The formal differential operator ( |2.5| ) is of the Sturm-Liouville type (cf. [0), and 
therefore we are allowed to apply the Sturm-Liouville theory to our problem. 

Mathematically, all the information about the differential operator h that is provided 
by the Sturm-Liouville theory (resolvent, spectrum, spectral representation,...) is obtained 
from the generalized eigenvalue equation 

hx{r-E)=(^-^^ + V{r)^x{r-^E) = Ex{r-E), EeC, (2.6) 

subject to different boundary conditions. From a physical point of view, Eq. ( p.6|) is the time- 
independent Schrodinger equation. As mentioned in the introduction, the monoenergetic 
eigensolutions of ( p.6| ) are not in general square integrable, i.e., they are not in the Hilbert 
space. Those monoenergetic eigensolutions will be associated to antilinear functionals Fe G 
$x by 



FE{ip)= / drip{r)x{r;E). (2.7) 
Jo 

These functionals are generalized eigenvectors of the Hamiltonian H, 

H''Fe = EFe, (2.8) 

or more precisely, 

{iflH^'lFE) = {Hip\FE) = E{ip\FE) , V</. e $ . (2.9) 



2.2. Self-Adjoint Extension 

Our first objective will be to define a linear operator on a Hilbert space corresponding 
to the formal differential operator h and investigate its self-adjoint extensions. Among all 
the possibilities, we shall choose the self-adjoint extension that fits spherically symmetric 
potentials. Later sections will deal with the spectral properties of this self-adjoint extension 
and with the RHS induced by it. 

The Hilbert space that is in the RHS of the square barrier potential is realized by the 
space -^^^([0, oo), dr) of square integrable functions /(r) defined on the interval [0, oo). In this 
section, we find a subdomain 'D{H) of this Hilbert space on which the differential operator 
h is self-adjoint. This domain must be a proper dense linear subspace of -Z^^([0, oo), dr). The 
action of h must be well-defined on T>{H), and this action must remain in L'^{[0, oo), dr). We 
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need also a boundary condition that assures the self-adjointness of the Hamiltonian. Among 
all the possible boundary conditions that provide a self-adjoint extension (see Appendix 
we choose /(O) = 0. These requirements can be written as 



f{r)eL\[0,oo),dr), 
hfir)eL'{[0,oo),dr), 
f{r)eAC'[0,oo), 
/(0) = 0, 



(2.10a) 
(2.10b) 
(2.10c) 
(2.10d) 



where AC'^[0,oo) denotes the space of functions whose derivative is absolutely continuous 
(see Appendix 0). Condition (|2.10a|) just means that the wave functions are square nor- 
malizable. Condition ( p.lObp assures that the action of h on any /(r) G T>{H) is square 
integrable. Condition ( |2.10c|) is the weakest condition sufficient for the second derivative 
of /(r) to be well-defined. In our example, this condition implies that /(r) and /'(r) are 
continuous at r = a and at r = 6. Condition ( p.lOd|) selects the self-adjoint extension needed 
in physics. 

The reason why we choose ( |2.10d| ) is the following: in physics ||l|,p^-p^, the set of 
boundary conditions imposed on the Schrodinger equation ( |2.(j| ) always includes 

(2.11a) 
(2.11b) 



xiO;E) = 0, 

x{r; E), and x'{^] E) are continuous at r = a and at r = 6 . 



Condition ( |2.11b| ) is implied by ( ^.lOc ), so we just need to recover ( 2.11a| ). This is why we 
impose (|2.10d| ). 

The set of conditions ( |2.10| ) leads to the domain 



V{H) = {/(r) I f{r)M{r) G L\%^),dr), f{r) G ^^^[0,00), /(O) = 0} . 



(2.12) 



In choosing ( p. 121) as the domain of our formal differential operator /i, we define a linear 
operator H by 



Hf{r) := hf{r) 



2m dr'^ 



+ V{r) f{r) , /(r) G V{H) . 



(2.13) 



2.3. Resolvent and Green Functions 

The Green function is the kernel of integration needed to write the resolvent of H as an 
integral operator, 

{E - H)-' f{r) = / G(r, s; E)f{s) ds . (2.14) 
Jo 

The procedure to compute the Green function of our operator ( p.l3| ) is explained in |21[] (see 
also l2^). For the sake of completeness, we include in Appendix the theorem that is used 
to calculate G(r, s; E). 
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The expression of the Green function will be given in terms of eigenfunctions of the differ- 
ential operator h subject to different boundary conditions (see Theorem 1 of Appendix 0). 
We shall consider three regions of the complex plane and compute the Green function for 
each region separately. In all our calculations, we will use the following branch of the square 
root function: 

-.{E eC\ - TT < arg(E) < tt} I — >{E eC\ -n/2< arg(E) < n/2} . (2.15) 



Region ^{E) < 0, Q{E) ^ 
For ^{E) < 0, '^{E) 7^ 0, the Green function (see Theorem 1 of Appendix 0) is given by 



2m/h? xir;E)0{s;E) 



G{r,s;E) 



^~2m/h? E 2Jz{E) 

2m I r? x(s;E)e(T;E) 
y/-2m/h2 E 2J3{E) 



r < s 
r > s 



The eigenfunction x(r; E) satisfies the Schrodinger equation (| 
tions 

x{0;E) = o, 

x{r-E)eAC\[0,oo)), 
x{r; E) is square integrable at . 

The boundary conditions ( |2.17| ) can be written as 



lk{E) < 0, %E) ^0. (2.16) 
and the boundary condi- 





X(0; 


E) 


= 0, 


x(« 


-0; 


E) 


= x(« + 0;^), 


X'(a 


-0; 


E) 


= x'(a + 0;^), 




-0; 


E) 


= x(& + 0;i?), 


X'{b 


-0; 


E) 


= x{h + Q;E), 




x{r; 


E) 


is square integrable at , 



(2.17a) 
(2.17b) 
(2.17c) 



(2.18a) 
(2.18b) 
(2.18c) 
(2.18d) 
(2.18e) 
(2.18f) 



and lead to 



x(r;E) 



— e 



-^Er 



Ji(E)ev/^^-^»)^ + J,(E)eV-^^^-^°^- 
J,{E)e^-"^'''- + ME)e-^' " 



'-Er 



< r < a 
a < r < b 
b < r < oo 



(2.19) 



The functions J'l — J'4 are such that xij\E) satisfies the boundary conditions (|2.18|) , and 
their expressions are given by Eq. ( |B4| ) of Appendix 

The eigenfunction B(r;i?) satisfies the Schrodinger equation (|2.6| ) and the boundary 
conditions 



e(r;E) e AC2([0,oo)), 

Q{r]E) is square integrable at oc 



(2.20a) 
(2.20b) 
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The boundary conditions ( |2.20| ) can be written as 



and lead to 



e{a-0;E) = e{a + 0;E), 
e'{a-0;E) = e'{a + 0;E), 

6(6 -0;E) = e{b + 0;E), 
e'{b-0;E) = e'(6 + 0;E), 

Q{r; E) is square integrable at oo 



Ai{E)ey''^'''' + A2{E)e 



Q{r]E) = I A^{E)e\ 



Ai{E)e 



'-{E-Vo)r 



<r < a 
a < r < b 
b < r < oo 



(2.21a) 
(2.21b) 
(2.21c) 
(2.21d) 
(2.21e) 



(2.22) 



The functions Ai — A^ are such that 6(r; E) satisfies the boundary conditions ( p.21j ), and 
their expressions are given by Eq. (|BB|) of Appendix O. 



Region ^{E) > 0, Q{E) > 
When ^{E) > 0, ^{E) > 0, the expression of the Green function is 

'im./n^ x{r;E)e+{s;E) ^ ^ ^ 



y/2m/n? E "iiME) 



(2.23) 



The eigenfunction x(r; E) satisfies the Schrodinger equation (2^) and the boundary condi- 
tions (PT7| ), 



x{r]E) 



sm 



^Er) 



J3(i^)eV^^^ + J4(E)e- 



< r < a 

a < r <b 
6 < r < oo 



(2.24) 



The functions Ji — are determined by the boundary conditions (|2.18| ), and their expres- 
sions are hsted in Eq. (p8| ) of Appendix 0. 

The eigenfunction B+(r;£') satisfies the Schrodinger equation (2^) and the boundary 
conditions ( p.2q ), 



At{E)e 



'-Er 



+ At{E)e 



^Et 



<r < a 



0+(r; E) = \ ^+(E)eV1^(^-^°)'- + ^|(E)e-V^(^-^°)^ a<r<b 

b < r < oo . 



(2.25) 



The functions Ai — A'^ are determined by the boundary conditions (|2.21|) , and their ex- 
pressions are hsted in Eq. (p9|) of Appendix |B[ 
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Region ^{E) > 0, Q{E) < 
In the region ^{E) > 0, '^{E) < tlie Green function reads 

{_ 2m/^ xir;E)e-is;E) 
_\^xtoaej^,.^^ S(i5)>0, O(£)<0. (2.26) 

The eigenfunction x(r; ii^) is given by ( ^■24D . The eigenfunction 0_(r;i?) satisfies the 
Schrodinger equation ( |2.6| ) and the boundary conditions (|2.20| ), 



^r(-g)e'v ^^' + >4 ^(^)e~V^^' 0<r<a 

0-(r; i?) = <! ^3 (E)eV¥(^'- + ^j(E)e"V¥(^'- a < r < 6 (2-27) 

i y < r < oo . 

The functions Ai — are such that B_(r; i?) and its derivative are continuous at r = a 
and at r = 6. Their expressions are hsted in Eq. (pil|) of Appendix |B[ 

2.4. Diagonalization of H and Eigenfunction Expansion 

In the present section, we diagonalize the Hamiltonian and construct the eigenfunction 
expansion generated by the eigenfunctions of the differential operator h. In order to do so, we 
compute the spectrum of H and then construct a unitary operator U that transforms from 
the position representation into the energy representation. We will see that the spectrum of 
H is the positive real line [0, oo). In the energy representation, H will act as the multipli- 
cation operator, the Hilbert space will be realized by L^{[0, oo), dE) and the domain of the 
Hamiltonian will be realized by the maximal domain on which the multiplication operator is 
well-defined. On our way, we shall take advantage of some theorems of the Sturm-Liouville 
theory that are proved in Ref . . For the sake of completeness, we include those theorems 
in Appendix |C[ 

24.1. Spectrum of H 

We first compute the spectrum Sp(if) of the operator H by applying Theorem 4 of 
Appendix y (see also [^). Since H is self-adjoint, its spectrum is real. The spectrum is the 



subset of the real line on which the Green function fails to be analytic. This non-analyticity 
of G(r, s; E) will be built into the functions 0^{E) that appear in Theorem 4 of Appendix 0. 

From the expression of the Green function computed above, it is clear that the subsets 
(—00, 0) and (0, 00) should be studied separately. We will denote either of these subsets by 
A. 

Subset A = (-00,0) 

We first take A from Theorem 4 of Appendix |C| to be (—00, 0). We choose a basis for the 
space of solutions of the equation ha = Ea that is continuous on (0, 00) x A and analytically 
dependent on E as 
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Bi{E)ey ~^^" + B 2{E)e 



I 2rr 



Er 



iEr 



<r < a 
a < r < b 
b < r < CO 



a2{r;E) = Q{r;E). 



(2.28a) 



;2.28b) 



The functions Bi — B4 are such that o"i(r; E) and its derivative are continuous at r = a and 
at r = 6. Their expressions are hsted in Eq. ( |C9|) of Appendix |C[ The function Q{r] E) is 
given by Eq. (g]^. 
Obviously, 



X(r; E) = JsiEMr; E) + ME)a2{r-, E) , 
which along with Eq. (p.l6| ) leads to 



G{t,s-E) 



2m /n^ 1 
-2m/f9E^ 



a,{r-E) + B^a,{r-E) 



(r2{s;E) 



ME) 

r<s, ^{E) <0,'^{E) ^ 



Since 



a2{s;E) = (T2{s]E) , 



we can write Eq. (|2.3(]| ) as 



Gir,s;E) 



2m/ h"^ 



-2m/h^E 2 



, ME) 



ai(r; E)a2{s- E) + a2{r; E)a2{s; E) 

ME) 

r<s, ^{E) <0,'^{E) ^0. 



On the other hand, by Theorem 4 in Appendix y we have 

2 

G{r,s-E) = Y,0;,{E)a,{r-E)aj{s;E), 



r < s 



By comparing Eqs. (|2.32|) and ( p.33| ) we see that 

g 2m/ h"^ 1 



0-AE) 



Q 2m /rf 1 J4{E) 

y/-2m/h? E'^ ME) 



^{E) < 0, '^{E) ^ 



(2.29) 



(2.30) 



(2.31) 



(2.32) 



(2.33) 



(2.34) 



The functions 0^-{E) are analytic in a neighborhood of A = (—00, 0). Therefore, the interval 
(— cxD, 0) is in the resolvent set Re(if) of the operator H . 



Subset A = (0, 00) 
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Now we study the case A = (0, oo). In order to be able to apply Theorem 4 of Ap- 
pendix 0, we choose the following basis for the space of solutions of ha = Ea that is 
continuous on (0, oo) x A and analytically dependent on E: 



ai{r;E) = x(r; 



cos 



2m 



(J2{r;E) = { 



Ci{E)e 



Er) 



+ C2{E)e 



[Cs{E)e^^''' + C,{E)e 



iEr 



<r < a 

a < r < b 
b < r < oo 



(2.35a) 



(2.35b) 



The functions Ci — C4, whose expressions are given by Eq. ( |C10|) of Appendix y, are such 
that (72 and its derivative are continuous at r = a and at r = 6. The eigenfunction x(r; E) 
is given by Eq. ( |2.24| ). 

Eqs. (|;25D, (^:27D and (pSj) lead to 



E) = -^Mr; E) + §^a,{r; E) 



W{E) 



W{E) 



and to 



e.{r;E) = ^a,{r;E)-§^a,{r;E), 



W{E) 



W{E) 



where 



W{E) = J,{E)Cs{E) - J^{E)C,{E) . 
By substituting Eq. (|2]3|) into Eq. (|2]2|) we get to 

2m /h^ 1 



(2.36) 



(2.37) 



(2.38) 



G{t,s-E) 



a2{r;E) 



^2m/h?E 2iME) [ W{E) ' ' W{E) 

^{E) > 0,$5(E) > 0, r > s 



(2.39) 



By substituting Eq. (|237| ) into Eq. (|2]2|) we get to 

2m/n2 1 



Gir,s;E) 



^2m/h?E 2iJ^{E) 



C,{E) UE) , 
— ; — rO"i r: E) ; — rO"2 r; E) 

^{E) > 0,'^{E) < 0, r > s 



ai{s;E) 



Since 



Eq. (I^D leads to 



cTi(s; E) = ai{s] E) , 



(2.40) 



(2.41) 



Gir,s;E) 



2m/h^ 



^2m/h^E 2iJ^{E) 



C,{E) 
'W{E) 



ai{r- E)ai{s; E) + 



ME) 



W{E) 

^{E) > 0,$5(E) > 0, r > s 



a2{r; E)a,{s; E) 

(2.42) 
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and Eq. ( p.4q ) leads to 



Gir,s;E) 



W{E) 



a^{r-E)ai{s-E) 



J?XE) 



W{E) 

^{E) > 0, 55(E) < 0, r > s 



a2{r; E)ai{s; E) 

(2.43) 



The expression of the resolvent in terms of the basis ai, a2 can be written as (see Theorem 4 
in Appendix |C|) 



G(r, s; E) = ^ e+{E)a,{r- E)a,{s- E) , r > s , 



(2.44) 



By comparing (|2.44|) to (|2.42|) we get to 



2m/h? 1 -CjjE) 
y/2m/h^E^iME)W{E) 
2m /n? 1 1 



^{E) > 0, "^{E) > 0, 



(2.45) 



y/2m/h^ E 2i W(E) 

By comparing (|2.44|) to (|2.43|) we get to 



2m I 1 Cz(E) 

y/2m/h^ E 2i Jz{E)W{E) 
2m I 1 1 



3?(E) > 0, $5(E) < 0, 



(2.46) 



From Eqs. (|2.45|) and ( p.46|) we can see that the measures pi2, P21 and P22 in Theorem 4 of 
Appendix y are zero and that the measure pn is given by 



1 rE2—S 

pnm, E,)) = Jim hm — / [et,{E - re) - et,{E + ^e)] dE 

1 2mlh^ 1 



which leads to 



p{E)^p,,{E) 



E, 47r ^2m/h^E Jz{E)J^{E) 
1 2m /h^ 1 



dE. 



(2.47) 



47r ^2m/t9E \ME)\^ 



E e (0,cx)). 



[2.48) 



The function O^-^^E) has a branch cut along (0, cxd), and therefore (0, cxo) is included in 
Sp(if). Since Sp(if) is a closed set, Sp(if) = [0,oo). Thus the resolvent set of H is 
Re{H) = C - [0,00). 



2.4-2. Diagonalization and Eigenfunction Expansion 

We are now in a position to diagonalize the Hamiltonian. By Theorem 2 of Appendix 0, 
there is a unitary map U defined by 
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U : L2([0, oo), dr) i — > L2((0, oo), p{E)dE) 

/(r) ^ f{E) = {Uf){E) = / rfr/(r)^(^ (2.49) 

Jo 

that brings 'D{H) onto the space 

= {/(E) GL2((0,oo),p(E)rfE)| / dEE''\f{E)\^p{E)<cx^}. (2.50) 

Eqs. ( |2.49|) and ( p.50|) provide a p-diagonahzation of H. If we seek a 5-diagonahzation, i.e., 
if we seek eigenfunctions that are 5-normahzed, then the measure p{E) must be absorbed 
by the eigenfunctions and by the wave functions.^ This is why we define 

a{r-E):=^/^x{r;E), (2.51) 

which is the eigensolution of the differential operator h that is 5-normahzed, and 

f{E) := ^KmKm , J{E) e L\iO,oo),piE)dE) , (2.52) 

and construct the unitary operator 

U : L2((0, oo)), p{E)dE) i — > L'^{{0, oo), dE) 

7^ f{E) = Uf{E) := y^)f{E) . (2.53) 

The operator that (5-diagonahzes our Hamiltonian is U := UU, 

U : L\[0, oo)), dr) i — > L\{0, oo), dE) 

f^Uf:= f. (2.54) 

The action of U can be written as an integral operator, 

f{E) = {Uf){E)= drfir)'^i^, f{r)EL\[0,oo),dr). (2.55) 
Jo 

The image of T>{H) under the action of U is 

poo 

V{E) ■=UV{H) = {f{E) E L\{0,oo),dE)\ / E^\ f {E)\^dE < oo} . (2.56) 

Jo 

Therefore, we have constructed a unitary operator 

U : V{H) C L\[0, oo), dr) i — > V{E) C L2((0, oo), rfE) 

f^f = Uf (2.57) 

that transforms from the position representation into the energy representation. The op- 
erator U diagonalizes our Hamiltonian in the sense that E = UHU~^ is the multiplication 
operator. 



^The meaning of the (5-normalization of the eigenfunctions will be explained in Section 2.9. 



13 



E : V{E) C L\{0, oo), dE) i — > L\{0, oo), dE) 

f^{Ef){E):=Ef{E). (2.58) 

The inverse operator of U is given by (see Theorem 3 of Appendix y) 

POO 

f{r) = {U-'f){r)= dEf{E)a{r,E), f{E)eL\{0,^),dE). (2.59) 
Jo 

The operator U^^ transforms from the energy representation into the position representation. 
The expressions (|2.55|) and (|2.59|) provide the eigenfunction expansion of any square 



integrable function in terms of the eigensolutions a{r; E) of h. 

The unitary operator U can be looked at as a sort of generahzed Fourier transform: the 
Fourier transform connects the position and the momentum representations. U connects the 
position and the energy representations. The role played by the plane waves e"*^"^' (which 
are generalized eigenfunctions of the operator —id/dx) is here played by the cr(r; E) (which 
are generalized eigenfunctions of the differential operator h). Therefore a{r;E) = {r\E), 
which are the (5-normalized eigensolutions of the Schrodinger equation, can be viewed as 
"transition elements" between the r- and the i?-representations. 

The label / of the functions in the position representation is different to the label / of the 
functions in the energy representation, because they have different functional dependences. 
Similar considerations apply to the Hamiltonian, the domain, the resolvent, etc. This is not 
the standard practice in the physics literature, where different representations are usually 



identified and labeled by the same symbol (see, for instance, p2|- p^p!8| ). 



2.5. The Need of the Rigged Hilbert Space 

The Sturm-Liouville theory only provides a domain T>{H) on which the Hamiltonian H 
is self-adjoint and a unitary operator U that diagonalizes H. This unitary operator induces 
a direct integral decomposition of the Hilbert space (see 0,^), 



Hi — >un = n = ® / n{E)dE 

Jsp{H) 

f^Uf^{f{E)}, f{E)en{E), (2.60) 

where Ti. is realized by ( [0, oo), (ir), and H is realized by L"^ {[0, oo), dE). The Hilbert 
space Ti-lE) associated to each energy eigenvalue of Sp(if) is realized by the Hilbert space 
of complex numbers C. On Ti, the operator H acts as the multiplication operator, 

Hf^UHf = {EfiE)}, feV{H). (2.61) 
The scalar product on Ti. can be written as 



/ (f{E),g{E)] dE, (2.62) 
where the scalar product ( ■ , ■ )_£; on H^E) is the usual scalar product on C, 
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[f{E),g{E)^^ = f{E)g{E). (2.63) 

As we shall explain below, the direct integral decomposition does not accommodate 
some of the basic requirements needed in Quantum Mechanics. These requirements can be 
accommodated by the RHS. 

One of the most important principles of Quantum Mechanics is that the quantity {ip, Hip) 
should fit the experimental expectation value of the observable H in the state (f. However, 
{(p, Hip) is not defined for every element in Ti, but only for those square normalizable wave 
functions that are also in V{H). Therefore, not every square normalizable function can 
represent a "physical wave function", but only those that are (at least) in V{H). Another 
fundamental assumption of quantum physics is that the quantity 

disp^H=iip,H'ip)-iip,Hip)' (2.64) 

represents the dispersion of the observable H in the state ip, and that 



A^H = Jdisp^H (2.65) 



represents the uncertainty of the observable H in the state ip. The quantities ( 2.64|) and 



( ^.65| ) are not defined for every element of the Hilbert space either. Therefore, we would 
like to find a subdomain $ included in "D^H) on which the expectation values 

((^,if», n = 0,l,2,..., (^e* (2.66) 

are well-defined. 

Another important requirement of Quantum Mechanics is that algebraic operations such 
as the sum and multiplication of two operators are well-defined. In the HS formalism, 
these algebraic operations are not always well-defined because the domains on which these 
operators are self-adjoint do not remain stable under their actions in general. In fact, much 
of the trouble of the HS formalism comes from domain questions. In our case, the domain 
T>{H) in ( |2.12| ) does not remain stable under H. We therefore would like to find a subdomain 



$ included in T>[H) that remains stable under the action of H and all of its powers, 

— n = 0,1,2,... (2.67) 



One can see that if Eq. (|2.67| ) holds, then the expectation values ( |2.66| ) are well-defined for 



each ip in i.e., if the domain $ remains stable under the action of if, then the expectation 
values of H in any state ip & ^ are well-defined. 

In Quantum Mechanics, it is always assumed that for each E e Sp(ii) there is a Dirac 
ket \E) such that 

H''\E) = E\E) (2.68) 



and such that the Dirac basis vector expansion ( p..3| ) holds. Equation ( p.68| ) has no solution in 



the Hilbert space when E belongs to the continuous part of the spectrum of the Hamiltonian. 
In fact, Eq. ( p.68|) has to be related to the equation 



(x|/i^|E) = , (2.69) 
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which in the radial representation reads 

ha{r; E) = Ea{r; E) , (2.70) 

where h is the differential operator ( p.5| ) and a{r; E) is the delta- normalized eigenfunction 
(|3l|). Since a{r;E) = {r\E) lies outside {[0 , oo) , dr) , i.e., 

cir |a(r;E)p = oo, (2.71) 



the corresponding eigenket \E), which is defined by 

|E) : $ I >C 



if I — > {ip\E) := ip{r)a{r; E)dr , (2.72) 
Jo 

should also lie outside the Hilbert space. We shall show that \E) is an element of 
In summary, what our mathematical framework should provide us with is: 

1. a dense invariant domain # on which all the powers of H and all the expectation 
values ( p. 661 ) are well-defined. 



2. smooth enough wave functions so that Eq. ( p.68| ) holds in the sense 



{V\H-\E) = E{^\E), (2.73) 

3. any wave function can be expanded by a Dirac basis vector expansion. 

In the direct integral decomposition formalism, there is not enough room for either of these 
three requirements. This is why we introduce the RHS. 



2.6. Construction of the Rigged Hilbert Space 

The first step is to make all the powers of the Hamiltonian well-defined. In order to do 
so, we construct the maximal invariant subspace V of the operator H, 

oo 

P := Pi P(if") . (2.74) 

n=0 

The space T> is the largest subspace of 'D{H) that remains stable under the action of the 
Hamiltonian H and all of its powers. It is easy to check that 

V = {y^e L\[0,oo),dr) I /iXr) G L\[0,oo),dr), h''^{0) = 0, (^("^a) = ^^''\b) = 0, 

n = 0,l,2,...; y;(r)GC°°([0,oo))}. (2.75) 



The conditions (6) = in ( p.75| ) come from taking the discontinuities of the 



potential V{r) at r = a and at r = 6 into consideration (cf. |^). 
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The second step is to find a subspace $ on which the eigenkets \E) of H are well-defined 
as antilinear functionals. For each E G Sp{H), we associate a ket \E) to the generalized 
eigenfunction a{r; E) through 

|E) : $ I ^ C 



if\ — ^{ip\E):= ip{r)a{r;E)dr = {Uip){E) . (2.76) 
Jo 

As actual computations show, the ket \E) in (|2.76|) is a generalized eigenfunctional of H if 
$ is included in the maximal invariant subspace of H, 

^CV. (2.77) 

Due to the non-square integrability of the eigenfunction a{r; E), we need to impose further 
restrictions on the elements of V in order to make the eigenfunctional \E) in Eq. (p.76|) 
continuous, 

poo 

/ dr |(r + l)"(/i + l)'"(^(r)|2 < oo, n, m = 0, 1, 2, . . . (2.78) 
Jo 

The imposition of conditions (|2.78| ) upon the space V leads to the space of wave functions 
of the square barrier potential, 

POO 

^ = {ipeV\ / rfr |(r + l)"(/i + l)'"(^(r)|^ < cx), n, m = 0, 1, 2, . . .} . (2.79) 
Jo 

On we define the family of norms 



\Mn,m--= \ / dr |(r + l)"(/^ + l)-^(r)|^ n, m = 0, 1, 2, . . . (2.80) 



The quantities (|2.8CI|) fulfill the conditions to be a norm (cf. Proposition 1 of Appendix 0) 
and can be used to define a countably normed topology r$ on $ (see 0), 

Lfa-^^P iff \\^a-V\\n,m ^0, W, m = 0, 1, 2, . . . (2.81) 

a— >oo a— >oo 

One can see that the space $ is stable under the action of H and that H is r^-continuous 
(cf. Proposition 2 of Appendix 0). 

Once we have constructed the space we can construct its topological dual as the 
space of r^-continuous antilinear functionals on $ (see [^) and therewith the RHS of the 
square barrier potential (for I = 0) 

$ C L2([0,oo),c/r) C . (2.82) 

The ket \E) in Eq. ( |2.76| ) is a well-defined antilinear functional on $, i.e., \E) belongs 
to (cf. Proposition 3 of Appendix 0). The ket l-E) is a generalized eigenvector of the 
Hamiltonian H (cf. Proposition 3 of Appendix 0), 

H''\E) = E\E), (2.83) 



17 



i.e., 



{ip\H''\E) = {Hip\E) = E{ip\E) , Vy?G$. (2.84) 

On the space all the expectation values of the Hamiltonian and all the algebraic 
operations involving H are well-defined, and the generalized eigenvalue equation (|2.84|) holds. 
As we shall see in the next section, the functions if of ^ can be expanded by a Dirac basis 
vector expansion. 



2.7. Dirac Basis Vector Expansion 

We are now in a position to derive the Dirac basis vector expansion. This derivation 
consists of the restriction of the Weyl-Kodaira expansions ( p.55| ) and ( |2.59| ) to the space 
If we denote {r\(f) = f{r) and {E\r) = a{r; E), and if we define the action of the left ket 
{E\ on If G $ as {E\ip) := (p{E), then Eq. ( |2.55| ) becomes 

POO 

{E\ip)= dr{E\r){r\ip) , cp e ^ . (2.85) 
Jo 

If we denote {r\E) = cr(r; E), then Eq. (|2.59| ) becomes 



POO 

{r\ip)= dE{r\E){E\ip) , cp e ^ . (2.86) 
Jo 

This equation is the Dirac basis vector expansion of the square barrier potential. In fact, 
when we formally write (|1.3| ) in the position representation, we get to ( p. 861) . 

In Eq. ( p.86|) , the wave function {r\(f) is spanned in a "Fourier-type" expansion by the 
eigenfunctions {r\E). In this expansion, each eigenfunction {r\E) is weighted by {E\ip) = 
(p{E), which is the value of the wave function in the energy representation at the point 
E. Thus any function (fir) = {r\(p) of $ can be written as a linear superposition of the 
monoenergetic eigenfunctions a{r;E) = {r\E). 

Although the Weyl-Kodaira expansions ( p.55| ) and ( |2.59| ) are valid for every element of 



the Hilbert space, the Dirac basis vector expansions (|2.85|) and (|2.86|) are only valid for 
functions </9 G $ because only those functions fulfill both 



lf{E) = {ip\E) (2.87) 

and 

{^\H-\E) = (HiflE) = E{ip\E) . C 



Another way to rephrase the Dirac basis vector expansion is the Nuclear Spectral 
(Gelfand-Maurin) theorem. Instead of using the general statement of we prove this 
theorem using the Sturm-Liouville theory (see Proposition 4 of Appendix]^). The Nuclear 
Spectral Theorem allows us to write the scalar product of any two functions ip^ip of ^ in 
terms of the action of the kets \E) on ip,ip: 

POO 

iip,^)= dE{ip\E){E\^/j), V<^,^G$. (2.89) 
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It also allows us to write the matrix elements of the Hamiltonian and all of its powers 
between two elements ip,il) of ^ in terms of the action of the kets \E) on ip,ip: 

POO 

iip,H''^)= dEE''{ip\E){E\tP) , V(^,V^ e = 1,2,... (2.90) 
Jo 



2.8. Energy Representation of the RHS 

In this section, we construct the energy representation of the RHS. Since the unitary 
operator U transforms from the position representation into the energy representation, the 
action of U on the RHS provides the energy representation of the RHS. 

We have already shown that in the energy representation the Hamiltonian H acts as the 
multiplication operator E. The energy representation of the space $ is defined as 

#:=f/$. (2.91) 

It is very easy to see that $ is a linear subspace of L'^ ([0, oo), dE). In order to endow $ 
with a topology r^, we carry the topology on $ into 

:= Ur^ . (2.92) 

With this topology, the space $ is a linear topological space. If we denote the dual space of 
$ by then we have 

^7^$^ = (f/$)^ = $^ . (2.93) 

If we denote \E) = U^\E), then we can prove that \E) is the antilinear Schwartz delta 
functional, i.e., \E) is the antilinear functional that associates to each function ip the complex 
conjugate of its value at the point E (see Proposition 5 of Appendix]^), 

|E) : $ I >C 

{lp\E) ■.= W)- (2.94) 

Therefore, the Schwartz delta functional appears in the (spectral) energy representation 
of the RHS associated to the Hamiltonian. If we write the action of the Schwartz delta 
functional as an integral operator, then the Dirac 5-function appears as the kernel of that 
integral operator. 

It is very helpful to show the different realizations of the RHS through the following 
diagram: 

H] ip{r) $ C L^([0, cx)), dr) C \E) position repr. 

[U iU iU"" (2.95) 

E; (p{E) $ C L2([0, oo), (ii?) C \E) energy repr. 

On the top line of the diagram ( p.95| ), we have the RHS, the Hamiltonian, the wave functions 
and the Dirac kets in the position representation. On the bottom line, we have their energy 
representation counterparts. 
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2.9. Meaning of the ^-Normalization of the Eigenfunctions 



In this section, we show that the 5-normahzation of the eigenfunctions is related to the 
measure dfi{E) that is used to compute the scalar product of the wave functions in the 
energy representation, 

POO 

(V^,^)= / ^)^|J{E)d|I{E) . (2.96) 
Jo 

We will see that if the measure in (|2.96| ) is the Lebesgue measure dE, then the eigenfunc- 
tions are ^-normalized, and that if the measure is p{E)dE, then the eigenfunctions are 
p-normalized. 

For the sake of simplicity, in this section we label the wave functions in the position and 
in the energy representation with the same symbol. With this notation, Eq. ( |2.86| ) reads 



(p{r)= / dE ip{E)a{r; E) , (2.97a) 
Jo 

/>oo 

i){r)= dEi){E)a{r]E). (2.97b) 
Jo 

Since ip{r),%lj{r) G L^( [0, oo), cir), their scalar product is well-defined. 



{^,ip) = I dr ip{r)^lj{r) . (2.98) 
Plugging ( p. 97] ) into (p.98| ), we obtain 



pOC POC POO 

dE dE' ip{E)^{E') / dr a{r; E)a{r; E') . (2.99) 
Jo Jo Jo 



If we use the Lebesgue measure dE, then the scalar product (p.96|) can be written as 



(^,V^)= / dEcp{E)^lj{E). (2.100) 
Comparison of (|2.99|) and ( p.lOO| ) leads to 

POO 

/ dr a{r; E)a{r; E') = 5{E - E') , (2.101) 
Jo 

i.e., the eigenfunctions cr(r; E) are 5-normalized. 

We now consider the case in which the eigenfunctions are p-normalized. If we use the 
measure dfj,{E) = p{E)dE, then the scalar product of (p and ip is given by 



{ip,^)= ipp{E)^p{E)p{E)dE, (2.102) 



where ipp{E) := ip{E)/^/p{E) and ^p{E) := ^{E)/^p{E). If we define ap{r;E) :-- 
a{r]E)/^p{E), then Eq. (|2^) reads 
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/■oo 

^{r)= / i^p{E)ap{r-E)p{E)dE, (2.103a) 

/>oo 

^{r)= / V^p(E)crp(r;E)p(E)o?E. (2.103b) 
Plugging Eq. ( pl03D into (|2:98|) , we obtain 



/•OO noo poc 

{^,^)= dE dE'^jE)i,p{E')p{E)p{E') drap{r-E)ap{r-E'). (2.104) 
^0 Jo Jo 

Comparison of {^J0§ and (|2]T0|) leads to 



/•°° 1 

dr ap{r- E)a,{r; E') = ^'^(^ - E') , (2.105) 

i.e., the eigenfunctions ap{r\E) are p-normalized. 

3. CONCLUSION 

In this paper, we have constructed the Rigged Hilbert Space 

* C L2([0,cx)),rfr) C (3.1) 

of the square barrier Hamiltonian and its energy representation 

$ C L2([0,oo),rfE) C (3.2) 

The spectrum of the Hamiltonian H is the positive real semiaxis. For each value E of the 
spectrum of H, we have constructed a Dirac ket \E) that is a generalized eigenfunctional 
of H whose corresponding generalized eigenvalue is E. In the energy representation, \E) 
acts as the antilinear Schwartz delta functional. On the space all algebraic operations 
involving the Hamiltonian H are well-defined. In particular, the expectation values of the 
Hamiltonian in any element of $ are well-defined. Any element of $ can be expanded 
in terms of the eigenkets \E) by a Dirac basis vector expansion. The elements of $ are 
represented by well-behaved functions, in contrast to the elements of the Hilbert space, 
which are represented by sets of equivalent functions that can vary arbitrarily on any set of 
zero Lebesgue measure. Therefore, it seems natural to conclude that a physically acceptable 
wave function is not any element of the Hilbert space but rather an element of the subspace 

The monoenergetic eigensolutions cr(r; E) of the time-independent Schrodinger equation 
are not square integrable, and therefore they cannot represent an acceptable wave function. 
Those eigensolutions has been used to define the eigenkets \E) that expand the physical 
wave functions G $ in a Dirac basis vector expansion. The eigenkets \E) belong to the 
space The eigensolutions cr(r; E) = {r\E) have been also used to construct the unitary 
operator U that transforms from the position representation into the energy representation. 

In our quest for the RHS of the square barrier potential, we have found a systematic 
method to construct the RHS of a large class of spherically symmetric potentials: 
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1. Expression of the formal differential operator. 

2. Hilbert space 7i of square integrable functions on which the formal differential operator 
acts. 

3. A domain T>{H) of the Hilbert space on which the formal differential operator is self- 
adjoint. 

4. Green functions (resolvent) of that self-adjoint operator. 

5. Diagonalization of the self-adjoint operator, eigenfunction expansion of the elements of 
Ti in terms of the eigensolutions of the formal differential operator, and direct integral 
decomposition of H induced by the self-adjoint operator. 

6. Subspace $ of 'D{H) on which all the expectation values of H are well-defined and on 
which the Dirac kets act as antilinear functionals. 

7. Rigged Hilbert space $ C 7i C 
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APPENDIX A: SELF-ADJOINT EXTENSION 

In this appendix, we list the possible self-adjoint extensions associated to the differential 
operator h. We first need some definitions (cf. [pl|]). 

Definition 1 By AC^([0, oo)) we denote the space of all functions / which have a 
continuous derivative in [0, oo), and for which /' is not only continuous but also absolutely 
continuous over each compact subinterval of [0, oo). Thus /*^^^ exists almost everywhere, and 
is integrable over any compact subinterval of [0, oo). At /' is continuous from the right. 

The space y4C^([0, oo)) is the largest space of functions on which the differential operator 
h can be defined. 

Definition 2 We define the spaces 
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nlilO, oo)) := {/ E AC'iiO, oo)) \f,hfE L'{[0, oo), dr)} (Al) 
^^([0, oo)) := {/ G AC\[0, oo) I /, G ^^([0, oo), dr)} (A2) 
7iQ([0, oo)) := {/ G 7-^^([0, oo)) | / vanishes outside some compact subset of (0, oo)} . (A3) 

Using these spaces, we can define the necessary operators to calculate the self-adjoint 
extensions associated to h. 



Definition 3 If h is the formal differential operator ( |2.5| ), we define the operators Hq 
and Hi on ^^^([0, oo), dr) by the formulas 

V{Ho) = nl{[0,oc)), Hof:=hf, feV{Ho). (A4) 
V{Hi) = nl{[0,oo)), HJ:=hf, feViHi). (A5) 



The operators i^o and Hi are sometimes called the minimal and the maximal operators 
associated to the differential operator h, respectively. The domain V{Hi) is the largest 
domain of the Hilbert space L'^{[0,oo),dr) on which the action of the differential operator 
h can be defined and remains inside L'^{[0, oo), dr). Further, H^ = Hi. 

The self-adjoint extensions of Hq are given by the restrictions of the operator Hi to 
domains determined by the conditions (see page 1306) 



/(O) + a/'(0) = 0, -oo<a<oo. (A6) 

These boundary conditions lead to the domains 

V^{H) = {/ G V{Hi) I /(O) + a /'(O) = 0} , -oo < a < oo . (A7) 

On these domains, the formal differential operator h is self-adjoint. The boundary condition 
that fits spherically symmetric potentials is /(O) = 0, i.e., a = 0. This condition selects our 
domain (|2.12|) , 

V{H) = V^^oiH) = {/ G V{Hi) I /(O) = 0} . (A8) 



APPENDIX B: RESOLVENT AND GREEN FUNCTION 



The following theorem provides the procedure to compute the Green function of the 
Hamiltonian H (cf. Theorem XIII.3.16 of Ref. ^): 

Theorem 1 Let H be the self-adjoint operator ( |2.13|) derived from the real formal 
differential operator ( |2.5| ) by the imposition of the boundary condition ( p.lOdp . Let ^ 0. 
Then there is exactly one solution x(r; E) of ih — E)a = square-integrable at and 
satisfying the boundary condition ( |2.10d[ ), and exactly one solution 9(r; E) of {h — E)a = 
square-integrable at infinity. The resolvent {E — H)^^ is an integral operator whose kernel 
G{r, s; E) is given by 
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G{r,s;E) 



2m x{r\E)@{s;E) 

W w(x,e) 

2m xis;E)e(r;E) 



W w(x,e) 
where W{x, ©) is the Wronskian of x and B 



r < s 



r > s . 



If we define 



k : = 
Q:- 



2m 



2m 



{E - K,) 



then the functions J{E) of Eq. (|2.19|) are given by 
ME) 



2 




2 




2 











A; 



,fca 



ME) 

and the functions ^^(E) of Eq. ( ^12^ ) by 



-1 - ^ 1 e 



-fca 



s 

A; 



s 

A; 



ME) 



—e 
2 



Q 



-kb 



fcfe 



1 

-e 
2 



—ka 



1 fca 
2 



1 - S j e'^'^A^iE) + ( 1 + I 1 e"V"Q«^4(^) 



The expression for the Wronskian of x ^i^nd G_ is 

W{x.^-) = -'^kJ^{E). 

If we define 



24 



then the functions J{E) of Eq. ( |2.24| ) are given by 



J^[E) = ie-*Q" ( sin(A;a) + ^ cos(A;a^ 
J^(E) = ^e^Q'^ (^sm{ka) - cos{ka] 



ME) = K 
ME) = \e 



—ikb 



ikb 



1 + I] e'^^'ME) + (l - I ) e-'^'UE) 



1 - I) e*«^Ji(E) + (l + f ) e-'^^'UE) 



and the functions A^{E) of Eq. ( p.25| ) by 



1 ^ 



-ika 



At{E) 



Aka 



The Wronskian of x ^^id G_|. is 

W{x,Q+)=2ikME). 
The functions ^^(-E) of Eq. ( p.27| ) are given by 



1 

— ( 

2 



-ifca 



k 

Q 

1 + 1^ e'^^^A^iE) + (^1 - I ) e-'^'^^4 (E) 



(B7a) 
(B7b) 

(B8a) 
(B8b) 
(B8c) 
(BSd) 



(B9a) 
(B9b) 
(B9c) 
(B9d) 

(BIO) 

(Blla) 
(BUb) 
(Bllc) 
(Blld) 



The Wronskian of x and G_ is 



W{x,e^) = -2tkJs{E). 



(B12) 
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APPENDIX C: DIAGONALIZATION AND EIGENFUNCTION 

EXPANSION 



The theorem that provides the operator U that diagonahzes H is (cf. Theorem XIII. 5. 13 
of Ref. m) 



Theorem 2 (Weyl-Kodaira) Let h be the formally self-adjoint differential operator 
( ^■51 ) defined on the interval [0, oo). Let H be the self-adjoint operator ( |2.13D . Let A be an 
open interval of the real axis, and suppose that there is given a set {ai{r]E), a2{r]E)} 
of functions, defined and continuous on (0, oo) x A, such that for each fixed E in A, 
{cri(r; E), cr2(r; E)} forms a basis for the space of solutions of ha = Ea. Then there exists 
a positive 2x2 matrix measure {pij} defined on A, such that 



1. the limit 



{Uf)i{E) = lim lim 



Od- 



f{r)ai{r;E)dr 



(CI) 



exists in the topology of -L^(A, {pij}) for each / in ^^^([0, oo), dr) and defines an iso- 
metric isomorphism U of i?(A)L^([0, oo), dr) onto L'^{A, {pij})] 

2. for each Borel function G defined on the real line and vanishing outside A, 

UV{G{H)) = {[/,] G L\A, {p,,}) I [GU e L\A, {p,,})} (C2) 

and 

{UG{H)fUE) = GiE){UfUE), z = 1, 2, i? e A, / G V{GiH)) . (C3) 



The theorem that provides the inverse of the operator U is (cf. Theorem XIII. 5. 14 of 
Ref. H) 

Theorem 3 (Weyl-Kodaira) Let H, A, {pij}, etc., be as in Theorem 2. Let Eq and Ei 
be the end points of A. Then 

1. the inverse of the isometric isomorphism U of i?(A)L^([0, oo), dr) onto L'^{A, {pij}) is 
given by the formula 

{U-'F){r)= lim lim H [Y.F,{E)a,{r-E)p,,{dE)\ (C4) 
where F = [Fi, F2] G L'^{A, {pij}), the limit existing in the topology of L^{[0, 00), dr); 
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2. if G is a bounded Borel function vanishing outside a Borel set e whose closure is 
compact and contained in A, then G{H) has the representation 

POO 

G{H)f{r)= fis)KiH,r,s)ds, (C5) 
Jo 



where 



K{H,r,s) = Y. I G{E)a,{s;E)aj{r;E)p,,{dE). (C6) 



The spectral measures are provided by the following theorem (cf. Theorem XIII. 5. 18 of 
Ref. my. 



Theorem 4 (Titchmarsh-Kodaira) Let A be an open interval of the real axis and O 
be an open set in the complex plane containing A. Let {ai{r]E), o'2{r; E)} be a set of 
functions which form a basis for the solutions of the equation ha = Ea, E E O, and which 
are continuous on (0, oo) x O and analytically dependent on E for E in O. Suppose that 
the kernel G{r, s; E) for the resolvent {E — H)~^ has a representation 



Gir,s;E) 



j:h=iOUE)a.ir;E)a,{s;E) 



^^^^,9-^{E)a,{r;E)a,{s;E) 



r < s , 
r > s , 



(C7) 



for all E in Re(if) fl O, and that {pij} is a positive matrix measure on A associated with H 
as in Theorem 2. Then the functions 6^ are analytic in Re(if) fl O, and given any bounded 
open interval {Ei, E2) C A, we have for 1 < j < 2, 



Pii((^i,^2)) = lim^^olim, .^+0^. 

= lim^^o lim,^o+ 2^ lE^+s^^tAE 



It) 



e-^{E + ie)]dE 
e±{E + ie)]dE. 



(C8) 



The functions B{E) of Eq. ( p.28a| ) are given by 



[ 1 + i ) 



B,{E) 
B,{E) 
B2{E) 



Q. 



1 



-ka 



1 ka 



1 + S j e'^'^BsiE) + - 2 1 e-'^^B^iE) 



(C9a) 
(C9b) 
(C9c) 
(C9d) 
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The functions C{E) of Eq. (|2.35b| ) are given by 

Ci(E) = ^e~'«" (cos{ka) - ^sm{ka)] 
C^{E) = -e'^" (cos(A;a) + ^sin(A;a^ 



1 

2^ 



—ikb 



ikb 



(ClOa) 
(ClOb) 
(ClOc) 
(ClOd) 



APPENDIX D: CONSTRUCTION OF THE RHS 



Proposition 1 Tlie quantities 

II 'y' II n.m 



dr \{r + lY{h + l)'^ip{r)\\ (/? G $ , n, m = 0, 1, 2, 



(Dl) 



are norms. 



Proof It is very easy to show that the quantities ( [OTI ) fulfill the conditions to be a 
norm, 

(D2a) 
(D2b) 
(D2c) 
(D2d) 



IIV^ + ^||n,m < \m\n,m + 
||aV^||n,m = |a| llV^lln.m, 
II V-' II n,"! — ! 

If ||v^||„,rrt = 0, then = . 



The only condition that is somewhat difficult to prove is (P2d| ): if ||v5||n,m = 0, then 

(l + r)"(/i + l)Xr) = 0, (D3) 

which yields 

(/i + l)'"(^(r) = . (D4) 



If m = 0, then Eq. ( p4D implies (y9(r) =0. If m = 1, then Eq. ( P^ ) implies that —1 is an 
eigenvalue of H whose corresponding eigenvector is (p. Since —1 is not an eigenvalue of H, 
must be the zero vector. If m > 1, the proof is similar. 

Proposition 2 The space $ is stable under the action of iJ, and H is r^-continuous. 
Proof In order to see that H is r$-continuous, we just have to realize that 

||-f^</'IU,m = \\{H + I)ip - ifWn^^ 

n,m n,m 

= llV'lln.m+l + ll'/^l|n,m • (D5) 
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We now prove that $ is stable under the action of H . Let To say that G $ is 

equivalent to say that ip &V and that the norms Hv^Hn.m are finite for every n,m = 0,1,2, .. . 
Since Hip is also in V, and since the norms ^ are also finite (see Eq. ( |D5| )), the vector 

Hip is also in 



Proposition 3 The function 
|E) : $ I — ^ C 

ip^{ip\E) :-- 



ip{r)a{r; E)dr = {Uip){E) . 



(D6) 



is an antilinear functional on $ that is a generalized eigenvector of (the restriction to ^> of ) 
H. 



Proof From the definition (p6|) , it is pretty easy to see that \E) is an antilinear func- 
tional. In order to show that \E) is continuous, we define 



MiE) := sup \a{r;E)\ . 

re[0,oo) 



(D7) 



Since 



\{if\E)\ = \UipiE)\ 



dr ip{r)a{r; E) 



dr \ip{r)\\a(r; E)\ 

POO 

<M{E) I dr\^{r)\ 
Jo 



= M{E) 
<M{E) 
= M{E)( 



OO 

dr 



1 + r 



dr ■ 



1 



:i + r)\^{r)\ 

1/2 



dr 



-^(^)ll^lli,o, 



;i + r)2 
1 

'l + r)2 



dr 1(1 + r)ip{r)[ 



1/2 



1/2 



:d8) 



the functional \E) is continuous when $ is endowed with the topology. 

In order to prove that \E) is a generalized eigenvector of H, we make use of the conditions 
( p.75| ) and ( 2.78 ) satisfied by the elements of 



{^\H-\E) = {Hv\E) 



dr 



d^ 



2m dr"^ 



+ V{r) ip{r)a{r;E) 



2m 



dip{r) 
dr 



a{r;E) 



+ 



2m 



ip{r) 



da{r;Ey 
dr 
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+ J dripir)i- — — + Vir)\a{r;E) 
= E{cp\E). (D9) 
Similarly, one can also prove that 

{^\{H-r\E) = E-{if\E). (DIO) 

APPENDIX E: DIRAC BASIS VECTOR EXPANSION 

Proposition 4 (Nuclear Spectral Theorem) Let 

^ C L\[0,oo),dr) C (El) 

be the RHS of the square barrier Hamiltonian H such that $ remains invariant under H 
and if is a r#-continuous operator on Then, for each E in the spectrum of H there is a 
generalized eigenvector \E) such that 

H''\E) = E\E) (E2) 

and such that 

M)= I dE{^\E){E\i;), \f^,ijE^, (E3) 

JSx>(H) 



and 



((/?,ii»=/ dEE''{ip\E){E\ij) , V(/?,V^ e $,n = 1,2,... (E4) 

Jsp{H) 

Proof Let (f and ip be in Since ?7 is unitary, 

(^,7A) = (f/^,f/^) = (^,^). (E5) 

The wave functions (p and ■?/' are in particular elements of L'^{[0,oo),dE). Therefore their 
scalar product is well-defined. 



(^,^)= / dEip{E)^{E). (E6) 

JSp{H) 

Since (p and ip belong to the action of each eigenket \E) on them is well-defined. 



{^\E) = ^{E) , (E7a) 
(^IV^) = ^{E) . (E7b) 

Plugging Eq. (|E^) into Eq. O and Eq. (p|) into Eq. (|E|), we get to Eq. (p]). The proof 
of fEif) is similar: 



30 



= / dE^{E){E''ij){E) 

Jsp{H) 

= / dE E''(p{E)i){E) 

Jsp{H) 

= I dEE''{ip\E){E\i^) . (E8) 



APPENDIX F: ENERGY REPRESENTATION OF THE RHS 

Proposition 5 The energy representation of the eigenket \E) is the antilinear Schwartz 
delta functional \E). 

Proof Since 

{(p\U-\E)^{U-'(p\E) 



(p{r)a{r; E)dr 



= m , (Fl) 
the functional U^\E) = \E) is the antilinear Schwartz delta functional. 
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